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A METHOD OF SOLVING LINEAR DIFFERENTIAL EQUATIONS* 

Br P. A. Lambert 

The object of this paper is to present a new method of solving ordinary 
linear diffential equations, which may frequently be applied with advantage 
when the coefficients of the equation are polynomials in the independent 
variable. 

Let the given differential equation be 

The method of solution proposed consists of the following steps : 

(a) Break up the function f into two parts, one of which, f lf equated to 
zero gives a differential equation which may be readily solved, and introduce 
a parameter t as a factor of the second part so that the given equation, 
yi +/ 2 = 0, is replaced by 

(2) A + «/, = 0. 
(&) Assume that the series 

(3) y = A + Bt + CP + Dt 3 - . • , 

where A, B, O, D, • • • are undetermined functions of x, satisfies (2). Sub- 
stitute the expression (3) for y in (2) and determine these functions by 
solving the differential equations formed by equating to zero the coefficients 
of successive powers of t in this identity. 

(c) Substitute the values of A, B, C, D, • • ■ in (3), and replace t by 
unity, and see whether 

* Read before the American Mathematical Society, October 30, 1909. 
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(4) y=A+B+ C + D... 

satisfies equation (1). 

It is to be noted that A is a solution of the equation /i = 0, and may con- 
sequently be frequently determined by a suitable choice of f v to contain n 
arbitrary constants so leading to the general solution of (1). If this is 
possible it is necessary only to obtain particular solutions of the equations 
for B, C,D, ■ ■ ■ 

Since the process above described is purely formal it is evidently neces- 
sary to see whether or not the series (4) actually satisfies (1) if that series 
contains an infinite number of terms, or if any of the functions A, B, C, • ■ • 
is given by an infinite series. 

The most advantageous method of breaking up the given equation into 
two parts must be determined by trial. However, if no term is separated 
into two parts, the number of possible methods of choosing/i and f± is never 
greater than 2 n , and if n is not large the best method may be selected without 
much difficulty. 

In the process of solving the differential equations which determine 
A, B, O, D, • • •, independent arbitrary constants are introduced until the 
number of arbitrary constants equals the order of the given differential 
equation. 

If the independent arbitrary constants are C lf C 2 , O s , • • • C„ the terms 
of the series (4) may be grouped so that (4) takes the form 

(5) y = C, y x + C 2 y 2 + C\ y s + • • • + C n y„ + T, 

where y x , y it y a , • • • y„ are independent solutions of the corresponding homo- 
geneous differential equation, and Fis a particular iutegral of (1). 

If another manner of breaking up the given differential equation makes 
series (4) a solution of (1), there may be determined a different set of in- 
dependent integrals and a different particular integral. 

If the general solution (5) of the differential equation contains infinite 
series, the limits of the regions of convergency must be determined by the 
usual methods. 

This method of solving differential equations is the result of an attempt 
to extend to differential equations the method employed by the author in the 
papers "On the solution of algebraic equations in infinite series." t 

t Bulletin American Mathematical Society, ser. 2, vol. 14, 1908, pp. 4»i7-477. 
Proceedings American Philosophical Society, vol. 47, 1908, pp. 111-134. 
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The incentive to make this attempt came from the following statement in 
an extract of a letter from Cauchy to Coriolis of January 29, 1837, published 
in the Comptes Iiendus of the Paris Academy. 

"Ainsi etendus, ces methodes s'appliquent avec un succes remarquable a 
presque tous les grands problemes d'analyse, a la resolution generale des 
equations, a l'integration des equations differentielles, a la mecanique 
celeste, etc." 

Cauchy describes the method applied to algebraic equations as follows : 

"Pour resoudre une equation partagcz son premier membre en deux poly- 
nomes d'une maniere quelconque, et supposez l'un de ces polynomes multiplie 
par un parametre que vous reduisez plus tard a l'unite." 

By this method the solutions of Bessel's equations, of Legendre's equa- 
tion, and of the differential equation of the hypergeometric series may be ad- 
vantageously determined. 

The method will be examplified by applying it to two differential 
equations. 



Example 1. Solve 



as? +«*v=i + *. 



Writing this equation in the form 



and assuming that 



there results 



g - (1 + x) + axhjt = 0, 



y = A + Bt + Ci?+ Dt 3 + 



cPA 
dx* 

(1 + x) 



+ 


&B 
dx 1 


+ 


ax 2 A 



t + 



d*C 



dx* 
+ axW 



< 2 + 



d*D 



dx* 
+ ax*C 



t 3 + ■ 



Equating to zero the coefficients of the successive powers of t in this 
identity, we have 
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CPA , 

M = 1+X ' 

cPC 



dx 2 



+ axW = 0, 



From this series of differential equations 
A= K t + K 2 x + ^ + ^, 



B = -R\ £L- K, 



ax 5 ax 6 ax" 1 



G = Zi ^nr~ a + K* 



3-4 2 4-5 1-2-5-6 2-3-6-7 ' 
a 2 x 8 „ a?x 9 a 2 x 10 



3.4.7-8^ 2 4.5-8-9 T 1-2.5-6-9-10 T 2-3-6-7.10.il' 

The law of formation of the successive coefficients .A, 5, C, D, • • • 
is evident, and the value of y, when I is made unity, becomes 

_ / ax* a*x 8 a 3 x 12 \ 

y = K x ^1 - — + g^TyTg - 3.4.7.U.H.12 + ■ • • J 

/ ax* a 2 x 8 a 8 x 12 \ 

+ A 2 x ^l - _ + 45g9 - 4. 5 . 8 .9. 12 .13 + " ' ' ) 

x 2 / ax* aht? \ 

+ "2 I 5T6 + 5-6.9 10 j 

a;" / ace* a 2 x 8 \ 

+ "6 V 6^7 + 6-7. 10-11 ~ "')' 

This value of y is composed of four infinite series, each convergent 
for all values of x, and is the complete solution of the given differential 
equation. 

Example 2. Solve 

*sg + (x + 2x')g-4y = 0- 
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Writing this equation in the form 



* 8 S+*S-% + ^$ = o, 



dx* 



dx 



dx 



and assuming that 



y = A + Bt + CP + DP + Efi + . 



there results 



x 2 — 
dx 2 

dA 

+ X dx~ 

- ±A 



+ x 2 



+ x 



+ 2x 2 



6?B 
dx 2 

dB 
dx 

4B 

dA 
dx 



t + X 2 


d 2 C 

dx' 2 


+ X 


dC 
dx 


— 


W 


+ 2x 2 


dB 

dx 



t 2 + • • • = 0. 



Equating to zero the coefficients of the successive powers of t in this 
identity, we have 

2 d 2 A dA 



dx 2 dx 

d?B dB . _ „ 
7, + x AB + 2x 

& tin* 



dA 



0, 



n d 2 C dO tn _ -dB n 

aar ax ax 



The solutions of the first equation for .A are .4 = A^x 2 and A = Iux~ 2 . 
(a) Substituting A = A'ix 2 in the next equation, 



d 2 £ 



d£ 



~ dx 2+ *dx-- 4Z?+4ir ^ = - 

A particular solution of this differential equation must be determined. 
The particular solution may be found by multiplying this equation by x 
and by two successive direct integrations. 



190 



LAMBERT 



[July 



The particular solution may also be found by the method of this paper 
as follows : Writing the equation in the form 



d"-B 

Tte 2 " 



+ x z±\ + 4/i> 3 - 4Bt = 
dx 



and assuming that 



B = B X + B.,t + B 3 l 2 + B^ + 



there results 






, dW } 


+ X dx* 


/It- 2 


X d£> 


M dx* 


dB, 

dx 


dB 2 

dx 


dB % 
dx 


+ 4A" 1 x 3 


- iB, 


- 4B. 2 



I* + . ■ . = 0. 



Equating to zero the coefficients of the successive powers of t in this 
identity, we have 

. dW t dB x . j- .. . 

x 2 — =-i + sb -j-i- + 4/Tia; 3 = 0, 
ax 2 ax 

d 2 i? 2 dB 2 .„ n 



x : 



X 



dx 2 dx 

, 2 ,PB, , _ di? 3 



dce ; 



^ + x 



dx 



4# 2 = 0, 



From this series of differential equations, we find the particular solutions 
B x = - * K^, B 2 = - (1) V 3 , A = - (J) *,■■■■ 



Hence 



* = - -J A > 3 -(£)** -G) 8 *.*-■■■=-£ *.* 



Substituting this value of B, 



.d*C dO .„ 2-3.4 
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A particular solution of this equation is 



- 5^ l ■ 



In like manner 






jv= (- i) n 



(n + 1) 2" 



5.6.7-8 • • • (n + 4) 
Substituting the values of A, B, C, D, ■ • •, and making < unity, 



+ C— IV' - 

+ ^ ' 5-G-7-8 • • • (n+ 4) 

which is convergent for all values of x and a solution of the given differ- 
ential equation. 

(b) Substituting A = A" 2 x~-, 

we have 

„ d-B dB „ . r , _ , - 

* 5i? + * ST ~ ** ~ *** = °- 

A particular solution of this equation, found by the same methods, i s 

B = -^K&-\ 

Substituting this value of B 

.(PC dC . n , 8 

A particular solution of this equation is 

It is evident that D = 0, E — • • • are particular solutions of the re- 
maining differential equations of the series. 
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Hence 



„=*(«-• -*«-»+*) 



is a solution of the given differential equation. 

The general solution of the given differential equation is the sum of 
these two independent integrals. 



Lehigh University, 

South Bethlehem, Pa. 
February, 1910. 



